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The FermifLAI collaboration recently reported the detection of starburt galaxies in the high en- 
ergy y-ray domain, as well as radio-loud narrow-line Seyfert 1 objects, which were previously un- 
usual extragalactic suspects for high energy emission. Motivated by the presence of sources close 
to the location of composite starburst/Seyfert 2 galaxies in the first year Fermi/LAT catalogue, 
we studied high energy y-ray emission from such objects, aiming at disentangling the emission 
of starburst and Seyfert activity. We analysed 1.6 years of Fermi/LAT data from NGC 1068 and 
NGC4945, which count among the brightest Seyfert 2 galaxies. We found an excess of high 
energy y-rays of 8.3(7 and 9.2a for 1FGL J0242.7+0007 and 1FGL J1305.4-4928, which are 
found to be consistent with the position of the Seyfert 2 galaxies NGC 1068 and NGC 4945, re- 
spectively. For both sources, we detect no significant variability nor any indication of a curvature 
of the spectrum. While the high energy flux of NGC 4945 is consistent with starburst activity, 
that of NGC 1068 is an order of magnitude above expectations, suggesting dominant emission 
from the active nucleus. We show that a leptonic scenario can account for the multi-wavelength 
spectral energy distribution of NGC 1068. High energy y-ray emission is thus revealed for the 
first time in a Seyfert 2 galaxy. If this result is confirmed in other objects, new perspectives would 
be opened up into the GeV band, with the discovery of a new class of high energy y-ray emitters. 



25th Texas Symposium on Relativistic Astrophysics - TEXAS 2010 
December 06-10, 2010 
Heidelberg, Germany 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Lieence. 



http://pos.sissa.it/ 



Seyfert 2 galaxies in the GeV band: jets and starburst 



J. -P. Lenain 



1. Introduction 



The Fermi/LAT collaboration recently reported the discovery of four radio-loud narrow-line 
Seyfert 1 galaxies [jlj, gj. Previously undetected at high energies, these objects thus constitute a new 
class of high energy emitters. Extending this idea to Seyfert 2 galaxies, we searched for high energy 
emission from two active galactic nuclei (AGN) among the closest and brightest in the X-ray sky, 
NGC 1068 and NGC4945. 

NGC 1068 is an archetypal Seyfert 2 galaxy, located at z = 0.003786, i.e. 14.4 Mpc away, 
harboring a hidden Seyfert 1 core. Antonucci & Miller proposed the AGN unification based on 
studies of this source [|J]. Both AGN and starburst activities are present in the central region of 
NGC 1068 [Q, |5j, the latter dominating the infrared emission of the broadband spectral energy 
distribution (SED) [§]. 

The presence of a Fermi/LAT source, 1FGL J0242.7+0007, in the region of NGC 1068 in the 
11-months Fermi/LAT catalogue (1FGL) [|7p, without any proposed counterpart in radio nor in y- 
rays, motivated us to analyse 1.6 years of data from Fermi/LAT to better assess the origin of this 
/-ray emission. 

NGC 4945, another Seyfert 2 galaxy at z = 0.001908, also exhibits starburst activity within 
its core. Its emission was detected up to soft /-rays, as observed with INTEGRAL/SPl, and was 
found to be Compton thick [§]. This source was reported as a /-ray emitter in the 1FGL catalogue, 
although the authors did not conclude on the starburst or AGN origin of the emission. 



2. Data analysis 

We analysed ~1.6yr of Fermi/LAT data, spanning from August 4, 2008 to March 15, 2010, 
from a region of interest of 10° in radius around NGC 1068, using the publicly available Science 
Tools 1 , and we followed the unbinned likelihood analysis scheme [01. 



Assuming a power-law shape for the source spectrum, the Test Statistic (TS) [10] of the like- 
lihood analysis is 68.6, corresponding approximately to a 8.3a source detection in the 100 MeV- 
100 GeV range. The best-fit location of the source is a 72 ooo = 2 h 42 m 46 s , <5/ 2 ooo = 0°2'14" with 
an error circle radius of 6.1' (68% confidence level, CL), and is fully compatible with the po- 
sition reported in the 1FGL catalogue. Given the angular distance between the Fermi source 
1FGLJ0242.7+0007 and NGC 1068, and the optical extension of 6.5' for the host galaxy, we 
propose that this Fermi/LAT source is actually associated with the Seyfert 2 galaxy NGC 1068. 

All the sources reported in the Fermi/LAT 11-months catalogue [R] within a radius of 15° 
around NGC 1068 were included in the likelihood analysis, and modelled with power-law spectra. 
For NGC 1068, we obtain Fi 00 MeV-iooGeV =(8.6 ±2.3) x 10~ 12 ergcnT 2 s" 1 and F = 2.31 ±0.13. 
We also tried to fit the data with a broken power-law or a log parabola, but this did not improve the 
likelihood. No significant variability was found in the data. 

We followed the same procedure for the analysis of Fermi/LAT data in the region of NGC 4945, 
for the source 1FGL J1305.4-4928. The likelihood analysis on NGC 4945 results in a TS of 85.3 
in the lOOMeV-lOOGeV band, equivalent to a 9.2a detection. The nominal position of the /-ray 
emission, located at a/2000 = 13 h 05 m 33 s , 5/2000 = — 49°26'44" with an error circle radius of 3.2' 



ittp://fermi. gsfc.nasa.gov/ssc/data/analysis/software/ 
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Figure 1: Relationship between SN rate, total gas mass and y-ray luminosity of NGC 1068, NGC 4945, 
NGC 253, M 82, the LMC and the Milky Way. 

(68% CL), is only 1.6' away from the position of NGC 4945. Assuming a power-law shape on the 
source energy spectrum, a photon index of T = 2.31 ± 0.10 and a flux of FiooMeV-iooGeV =(1.6 ± 
0.3) x 10~ n ergcm~ 2 s _1 are found. As for NGC 1068, the use of a broken power-law or a log 
parabola did not improve the likelihood, and no significant variability was found in the data, which 
are statistically consistent with a constant for both sources. 

The ISGRI data were reduced using the INTEGRAL Offline Scientific Analysis software 2 ver- 



sion 9.0. We have also checked the Swift/BNT [11] spectra of the two sources. The 18 months 



Swift/BAT spectra have been extracted from the BAT archive [|12| served by the HEAVENS source 
results archive (Walter et al., in prep.). 

3. Discussion 

As a first test to know whether the high energy emission of these two objects orginates from the 
starburst or the AGN component, a significant detection of variability would rule out the starburst 
hypothesis. However, due to the lack of statistics, no conclusions can be drawn about the variability 
of NGC 1068 or NGC 4945 from their light curves. 

We also compare their y-ray luminosity with those of the famous starburst galaxies NGC 253 



and M 82, both detected at high energies with Fermi/LAT [jl3J and at very high energies with 



ittp://www.isdc.unige.ch/integral/ 
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Figure 2: Spectral energy distribution of NGC 1068, including the Fermi/LAT spectrum. The black and red 
points are archival data from the NED, the red ones denote data taken from the central region of NGC 1068. 
For clarity, we only show the INTEGRAL IBIS/ISGRI data in blue in the hard X-rays. The EIC model for 
the outflow is shown in blue, and the corresponding SSC emission is shown in thin red and magenta lines for 
first and second order components, respectively. The thick red line shows the sum of the different emission 
components from the large outflow. The EIC component from the accretion disc is shown in green. 



H.E.S.S. (Tj] and VERITAS @, respectively. Both have a luminosity in the 100MeV-5 GeV 
band of the order of «10 40 ergs _1 , as detected with Fermi [13]. Computing the luminosities of 
NGC 1068 and NGC 4945 in the same energy band, for comparison, we obtain 1.7 x 10 41 erg s _1 
and 2.0 x 10 40 ergs _1 , respectively. High energy emission from starburst galaxies is thought to 
originate from the interactions of cosmic rays produced by supernova; with the ambient interstellar 
medium, and the y-ray luminosity Ly should then scale with the product of the supernova rate /?sn 
and the total gas mass M gas . We thus compare 7? S n, M gas and L y of NGC 1068 and NGC 4945 to the 
ones of NGC 253, M 82, the Large Magellanic Cloud (LMC) and the Milky Way (see Figure [[]), as 
well as their infrared and radio luminosities (see [16] for more details). These objects are the only 
extragalactic sources which are not AGN known to emit high energy y-rays. The y-ray luminosity 
and the SN rate of NGC 4945 are fully consistent with those of NGC 253 and M 82, hence even 
though this object is a composite starburst/AGN, its high energy y-ray emission could be explained 
only in terms of starburst activity. 

However, in the case of NGC 1068, /?sn is comparable to those of M 82 and NGC 253, but its 
radio and y-ray luminosities are higher by a factor ~10. This would suggest that its high energy y- 
ray emission is more likely dominated by the central AGN activity. This is also strengthened by the 
fact that radio maps of NGC 1068 clearly show a structured jet, on parsec- and kiloparsec-scales, 
shaped by the outflow from the central AGN [|n[ |8|. On the contrary, the radio morphology of 
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NGC 4945 shows an extended emission consistent with the optical morphology tracing the edge-on 
galaxy, consistent with a starburst emission. 

Assuming that the y-ray emission of NGC 1068 is due to the central AGN, we apply an emis- 
sion model developed for AGN which jet is misaligned to the line of sight [19]. A large, mildly- 
relativistic zone of the wind-like outflow, at a few tenth of parsecs from the core, could emit high- 
energy y-rays through external inverse Compton process (EIC) [20]. At such distances, the infrared 
photon energy density is still high enough to ensure a significant EIC emission while being not too 
important to prevent high optical opacity from pair production. In accordance with the data, no 
significant short-term variability is expected from such a large emitting zone. 

In the model, a compact, dense plasma radiates through synchrotron self-Compton (SSC) and 
EIC processes. Second order SSC emission is also accounted for. Figure ^] shows the results of our 
model, with synchrotron emission responsible for the radio emission, while the FermifLAI data are 
interpreted as EIC emission with the infrared emission providing the seed photons, from a multi- 
temperature blackbody. We also show the contribution from the SSC process, which is negligible 
compared to the EIC emission at the highest energies (see [^] for more details). The hard X-ray 
spectrum is due to EIC radiation from another relativistic population of leptons, within hot plasma 
located in the vicinity of the accretion disc. It should be noted that imaging atmospheric Cerenkov 
telescopes have the ability to strongly constrain the highest energy part of the particle energy dis- 
tribution, and hence provide insights on the acceleration processes at work in this object. If the 
leptonic population of particles extend at higher energies, a significant signal should be detectable 
from NGC 1068 with H.E.S.S., VERITAS, MAGIC or the future CTA (Cerenkov Telescope Array) 
observatory [pT|]. 

Figure ||] shows the SED of NGC 4945 including archival data, as well as the Fermi/LAI and 
INTEGRAL spectra. We also show for comparison the SED of the starburst galaxy NGC 253. The 
two sources clearly have very similar broadband SEDs, strengthening the idea that the high energy 
y-ray emission from NGC 4945 is of starburst origin. NGC 253 being detected at very high energies 
with H.E.S.S. jn]], NGC 4945 is thus expected to be also detectable in this energy range. 



4. Conclusion 



High energy y-ray emission is revealed and associated to the starburst/Seyfert 2 objects NGC 
1068 and NGC 4945, using Fermi/LAI data. Both source spectra are consistent with a power- 
law, and no variability is found in the data. Compared to M 82 and NGC 253, which high-energy 
emission is dominated by starburst activity, we find a too high y-ray luminosity of NGC 1068 to 
be explained only by starburst activity. We thus propose a leptonic scenario to interpret the high 
energy emission in terms of external inverse Compton process from an outflowing relativistic wind 
launched by the central AGN. 

If high energy y-ray emission due to AGN activity is confirmed in other Seyfert 2 galaxies, 
this would mark the discovery of yet a new class of high-energy y-ray emitters. New data from 
Fermi/LAI in the coming years will be extremely helpful to search for such emission. 
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Figure 3: Spectral energy distribution of NGC4945, including the FermilhPU spectrum (black points). For 
clarity, we only show the INTEGRAL IBIS/ISGRI data in blue in the hard X-rays. The model for the EIC 
component from the accretion disc is shown in green. We show in red the data of NGC 253 as taken from 
the NED, with its Fermi/LAT spectrum [jjjj] as well as the H.E.S.S. flux measurement [Q, for comparison. 
The luminosity axis on the right is given for NGC 4945. 
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